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I. INTRODUCTION 



Hawking's 1974 paper [jT[ triggered great interest in both the existence of the radiation 
which bears his name and speculations as to what his result implied for the validity of quan- 
tum field theory in a black hole background. Subsequent calculations, which showed that 
the phenomenon was robust 0, ^, Q and supported the view that the radiation appeared to 
be completely thermal, only added to this interest. Combining these results with specula- 
tions as to what happens when the black hole evaporates led Hawking and others to argue 
that the behavior of this system must be inconsistent with the unitary time evolution of 
the underlying field theory, since it starts out in a pure state and evolves into a thermal 
ensemble at a temperature T = l/(87rGM). 

A related issue, the so-called information paradox, arose after Bekenstein P] argued that 
a black hole of mass M has an entropy proportional to its area. Here, the main question is 
"What happens to information which has already crossed the horizon when the black hole 
evaporates?" . 

Both questions led to suggestions that something goes wrong with quantum field theory, 
even for a large semi-classical black holes, and that a careful study of questions related to 
these aspects of black hole physics would point the way to the theory which must replace it. 
Intrigued by this idea and convinced that at least some of these questions would be easier 
to analyze if we could canonically quantize the field theory in the background of a large 
Schwarzschild black hole, we decided to do just that. 

In this paper we present a comprehensive treatment of the work first discussed in an 
earlier Letter , wherein we canonically quantized a massless scalar field in the background 
of a large Schwarzschild black hole and derived the familiar Hawking results. In addition to a 
more detailed treatment of the problem we show that our approach leads to a different picture 
of Bekenstein's derivation of the entropy of a black hole and a surprising, but self-consistent, 
resolution of the information paradox. These results suggest that despite expectations, 
studying the problem of Hawking radiation at the purely semi-classical level is unlikely, in 
itself, to produce new insights into the question of how a quantum theory of gravity should 
behave. 
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II. PRELIMINARIES 



Before going further it is useful to review why a Hamiltonian formulation of the problem 
of a massless scalar field in the presence of a black hole background appears problematic. 
Let us begin by considering the problem of a massless scalar field with Lagrange density 

C=^[g>^''d^cP{x)d,cP{x)] (1) 

in the background of a Schwarzschild black hole of mass M. In the usual Schwarzschild 
coordinates the metric g^^ takes the familiar form 

ds^ = _(1 _ ^) de + (1 - — dr^ + r^dn\ (2) 

where we have set Newton's constant, G, to one. 

As is well known, the apparent metric singularity at r = 2M is a coordinate artifact 
and, as such, does not pose a problem. The true issue for canonical quantization is that 
we need to define a family of spacelike slices which foliate the spacetime in order to define 
initial data and form the Hamiltonian. Inspection of Eq.(0) shows that surfaces of constant 
Schwarzschild time change from spacelike to timelike at the horizon (r = 2M ) and so they do 
not fulfill our requirements. As we show in the next Section, changing to Painleve coordinates 
both eliminates the coordinate singularity at r = 2M and allows us to analytically define 
a satisfactory family of spacelike slices. Since, however, the form of the metric in Painleve 
coordinates is not well suited to simple canonical quantization, we need to introduce yet 
another transformation, to Lemaitre coordinates, to facilitate the quantization procedure. 
The Hamiltonian constructed in this way explicitly depends upon Lemaitre time, which is a 
specific manifestation of the general theorem that "the Schwarzschild metric does not admit 
a global timelike Killing vector field". Fortunately, having an explicitly time dependent 
Hamiltonian is no barrier to unitary time evolution; in fact, this is always the case for the 
interaction representation. The crucial requirement is not that the Hamiltonian is time 
independent, but rather that there exists a one-parameter family of unitary operators U{X) 
which satisfy the equation 

'^ = -^H{X)U{X). (3) 

The lesson we learn from the fact that the Hamiltonian explicitly depends upon time is 
that we shouldn't be looking for static quantities such as the vacuum state of the theory 



but rather for steady state phenomena such as the Hawking radiation. In a sense, once we 
have observed that the Hamiltonian is time dependent there is no longer a puzzle as to why 
Hawking radiation can exist. What remains is to fill in the details and show how explicit 
calculations in this canonical framework lead to Hawking's results. We do this in the next 
few sections. Once the equivalence of our discussion to previous approaches is established 
we turn to a discussion of what our approach has to say about the question of black hole 
entropy and the information paradox. 



III. COORDINATE SYSTEMS 



Although we are ultimately interested in the surfaces defined by constant Painleve time, it 
is convenient to begin by introducing Kruskal coordinates. First, because these coordinates 
make it particularly easy to draw null-geodesics (they are simply lines parallel to either the 
X or Y axes shown in Fig.|I|); second, because they allow us to easily compare surfaces of 
fixed Schwarzschild time to surfaces of fixed Painleve time. 

It is convenient to introduce dimensionless versions of r and t by rescaling r — > 2Mr 
and t 2Mt. Using these variables we introduce the Kruskal coordinates X and Y by the 
equations: 

Ar=(r-l)e^ 
In these coordinates the Schwarzschild metric takes the form 

^ 2 32 e^-'UXdY 2j^2 /rN 

ds^ = h r^dQ^. (5) 

r 

Eq. tells us that fixed Schwarzschild r is a hyperbola in the X, y-plane, as shown in Fig.|l|, 
and that a surface of fixed Schwarzschild time corresponds to a straight line X = |y| e*^ 
(such lines are not shown in Fig.|I|). 

Painleve coordinates are derived from Schwarzschild coordinates by making an r- 
dependent shift in Schwarzschild time; i.e., 

v^- 1 

This equation makes it easy to compute surfaces of fixed Painleve time. These surfaces are 
the almost horizontal curves shown in Fig.^and they clearly foliate the spacetime. Note that 
while Schwarzschild t and Painleve A differ by a function of r, two events having the same r 



t = A - 2Vr - In 



(6) 



r=2M 



Y 



FIG. 1: Only two of usual four X-Y Kruskal quadrants are plotted and the figure is rotated by 
45° to emphasize we are studying the region from r = to r = co. Vertical and horizontal lines 
are null- geodesies. The nearly, but not quite, horizontal curves are globally spacelike surfaces of 
constant Painleve time. The hyperbolas are surfaces of constant Schwarzschild r. Lines of constant 
Schwarzschild time would be straight lines originating at X = Y = 0. 

are separated by equal intervals of Schwarzschild or Painleve time. In Painleve coordinates 
the Schwarzschild metric takes the form 



Although Painleve coordinates are useful for defining a family of spacelike surfaces they 
are not well suited for canonical quantization because of the cross term dX dr in the metric 
and the fact that lines of constant r (we assume the angular variables 6 and are held fixed) 
are not everywhere timelike. A better coordinate system can be obtained by considering a 
set of curves (A, t(A)) whose tangent vectors, (1, dr{X)/dX), are orthogonal to the surfaces of 
constant Painleve time. Substituting this requirement into Eq.(|^ for the metric in Painleve 
coordinates we arrive at an equation which, upon integration, gives: 



where rsch is the value of Schwarzschild r at which each curve passes through the surface 
defined by A = 0. 






(8) 
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FIG. 2: This is the same plot of surfaces of constant Painleve time overlaid with curves showing 
r(A, r/) for two different initial values of r]. Note that all such lines finally intersect the spacelike 
curve r = 

We could use this coordinate system but to avoid dealing with factors of r^/^ it is conve- 
nient to make one more change of variables. This leads to Lemaitre coordinates, which are 
related to Painleve A and r by 

/ o X 2/3 /o X 2/3 

r(A,r,e.)=(rf4'-|A) =(|(^-A)) ■ (9) 

In Lemaitre coordinates the metric takes the form 

rfs' = -rfA2 + -l— (ir/' + r(A,r/)W. (10) 

It is manifestly free of coordinate singularities at r = 1, has no cross terms in dX and dr and 
allows a completely straightforward canonical quantization procedure. Fig. 2 shows lines of 
constant rj overlayed on the surfaces of constant Painleve time. 

IV. CANONICAL QUANTIZATION 

The Lagrangian for a massless scalar field theory has the general form given in Eq.(0). 
Since the metric in Schwarzschild, Painleve and Lemaitre coordinates is rotationally in- 
variant we can study the scalar field theory for each angular momentum mode separately. 
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Furthermore, we are free to expand the field 0(A, 77, 6, 0) in spherical harmonics in 9 and 
and restrict attention to the L = mode, since this exhibits all of the interesting behavior. 
The Lemaitre coordinate form of the L = scalar field Lagrangian is 

C^-y/^l [(5a0o (A, r))f - r (9,0o (A, r))f] (11) 



2 

where the determinant is 

v^ = r^/^ = ^(^-A). (12) 

Following the usual rules for canonically quantizing such a theory we see that the mo- 
mentum conjugate to the field is 

^o(A,r7) = ^^^^a,0o(A,r;), (13) 

and the canonical Hamiltonian is 

The commutation relations for 0o and ttq are 

[ko{\v,),4>oM)]^ -i^V-v')- (15) 

As we already pointed out, there is a one parameter family of unitary operators, C/(A), 
which satisfy the equation 

^U{X)^-iH{\)U{Xl (16) 
whose solution is the path ordered exponential 

U{\) = V [e-' lo '^^^W^ . (17) 

Given these operators we define fields at later times as 

0o(A, 77) = U{X) Mv) U^X), 7ro(A, 77) = U{X) iroiv) C/^(A). (18) 

It follows from the canonical commutation relations that these operators satisfy Heisenberg 
equations of motion of the form 

dx [{v - X)dx<f>o] - d, [{rj - A) r d,<f>o] = 0. (19) 
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Clearly we have two options open to us. The first is to diagonalize the Hamiltonian 
H{0) and then explicitly construct U{X). This is at best cumbersome. The second, more 
tractable option for a system with a time-dependent Hamiltonian, is to solve the Heisenberg 
equations of motion and compute all physical quantities by evaluating expectation values of 
the interesting time dependent operators in a fixed initial state. We will adopt the second 
approach and use the vacuum state of the Hamiltonian H{0) as our initial state. It will 
be apparent from the computations which follow that except for transient effects, it would 
make no difference if we chose as our initial state any state whose energy differed from the 
H{0) vacuum state energy by any finite amount. Note that there is a subtlety associated 
with the computation of the time-dependent Hamiltonian at 77 = A or r{X,ri) = 0. We will 
return to this point after our discussion of the geometric optics solution of the Heisenberg 
equations and the derivation of Hawking radiation. This issue is important and relates to 
the way in which the theory deals with the information paradox. 

Before discussing the solution of the Heisenberg equations of motion let us point out that 
it is simple to find all of the eigenstates of H{0) because it is just a free field Hamiltonian 
in disguise. To see this we only need to change variables, back to Schwarzschild r, using 
rj = {2/3)r^^'^ and then rescale the fields by 

vro(r) = v^7ri(r), 0o(r) = (20) 

This converts Eq . ([I^) to 

H{0) = \ dr (^7ri(r)2 + r\dr ^)^) , (21) 

which is the Hamiltonian of the L = mode of a free massless field in fiat space. One 
constructs the eigenstates of this Hamiltonian in the usual way by expanding the fields in 
terms of annihilation and creation operators: 




and defining the vacuum state |0) to be the state that is annihilated by all of the a^'s. 

V. A LOOK AHEAD 

At this point we have to deal with two important issues. The first is to check how the 
formalism just presented leads to the usual Hawking results. The second is to see how this 
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Hamiltonian approach modifies the way we think about the questions of black hole entropy 
and the information paradox. 

Since one of the main arguments for the usual interpretation of Hawking radiation is that 
a black hole is a thermodynamic system, before discussing the case of the black hole it is 
worth spending some time reviewing the simpler problem of a constantly accelerating mirror 
in a 1 + 1 dimensional spacetime. We include this discussion to emphasize that a manifestly 
non-singular quantum mechanics problem with a time-dependent Hamiltonian and unitary 
time evolution can produce, what looks to some observers, like purely thermal outgoing 
radiation. Furthermore, although our discussion will cover well worn ground, the technique 
we use to derive this result differs somewhat from the approach presented in Ref. 0. In 
particular, we will focus on a simple configuration space solution of the field equations, 
relating all measurements back to the initial surface on which we quantized the theory (i.e., 
at time A = 0). This allows us to get a better handle on the physical assumptions being 
made when we discuss the initial state. 

In Section|VI| we set up the problem of the moving mirror and discuss the solution of 
the field equations in the same way we will discuss them for a field theory in a black hole 
background. One reason for studying this problem in detail is that later, after discussing 
the eternal black hole problem as it is usually formulated, we will reformulate it in the spirit 
of the moving mirror problem in order to justify our choice of the vacuum state on physical 
grounds. 

After setting up the moving mirror problem we discuss what an Unruh thermometer 
would measure at late times, showing that the process of adiabatically turning on and then 
turning off such a device leads to a measurement of what appears to be a non- vanishing 
temperature. Proceeding in the same vein we give a detailed derivation of the energy flux 
passing through a fixed position in space and show that it appears to be thermal. We do this 
to emphasize the fact that a perfectly unitary quantum system can exhibit some apparently 
peculiar behavior if the Hamiltonian is time dependent. Technically this computation is 
similar to the computation one has to do for the black hole case. The most important result, 
for both cases, is that although the energy density of the outgoing radiation is divergent and 
therefore ill defined, the flux computation is free of divergences and unique. 

These arguments are followed by Section where we attempt to match our approach to 
the Bogoliubov transformation technique, often discussed in the context of black hole physics. 
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Here we show that while one can estabhsh such a connection, matching the Bogohubov 
approach onto a Hamiltonian formahsm that exphcitly works with finite times is somewhat 
unnatural. 

Finally, wc conclude our excursion into theories in flat space with a variant of the moving 
mirror problem where there is both a moving and a fixed mirror. This problem is interesting 
because it exhibits a peculiar feature of the late time problem; namely, that at late times, 
the fields over most of space depend only upon degrees of freedom which, on the original 
surface of quantization, are localized within an exponentially small region surrounding a 
single point. This strange behavior is really just a refiection of the causal structure of the 
problem. It is a feature of the Schwarzschild problem as well and will play a role in our 
discussion of Bekenstein entropy. 

With our foray into non-gravitational physics behind us, we turn our attention to the 
case of the Schwarzschild black hole. A brief statement of what we mean by the geometric 
optics approximation to the field equations is followed by a recapitulation of the Unruh 
thermometer and energy fiux computations for the black hole. Once we have shown that 
our approach reproduces the well known results, we turn to things which can be better 
discussed in this framework. 

The first benefit which follows from discussing the black hole problem between surfaces 
separated by finite times is that we can discuss a variant of the problem of Hawking radiation 
which provides a rationale for picking a particular vacuum state. For this purpose we consider 
the problem of an infalling mirror, i.e., a version of the black hole problem in which, up to a 
time Ao the black hole is surrounded by a refiecting sphere of large, fixed radius i?o, which at 
time Ao falls into the black hole along a Lemaitre timeline. By a reflecting sphere we mean 
that we assume that the fleld always vanishes on and inside the surface. Since the field is in 
a region of vanishingly small gravitational field for an infinite time in the past, at least if Rq 
is chosen large enough, we can argue that from the point of view of a local observer in this 
region it is sensible to assume that the system starts out in the vacuum state. We then show 
how, in a manner very reminiscent of calculations done for a self-assembling black hole, the 
Hawking radiation forms as the mirror approaches the horizon. 

Another problem which becomes more approachable because we work between finite 
times is the so-called back reaction problem. The back reaction problem is equivalent to the 
observation that the eternal Schwarzschild background is not consistent with the addition 
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of the scalar field theory, since the energy momentum tensor we compute for the scalar field 
has a uniquely defined non- vanishing flux term and the Einstein tensor for the Schwarzschild 
solution vanishes. Since we can discuss this issue for large black holes and finite times, during 
which the mass of the hole does not change much, we argue it is possible to ask and solve 
the question of what a truly self-consistent problem would look like. Wc emphasize that this 
is possible because we avoid the issue of what happens at infinite times in the future when 
the evaporation process becomes rapid and runs to completion. 

To discuss the question of black hole entropy we study a variant of the black hole problem 
in which we place a large reflecting surface around the black hole, but now we assume that 
the field theory exists only inside this surface. Wc then show that to an outside observer a 
body constructed in this way appears to be a thermal system with the familiar Bckcnstein 
entropy, however when we look inside we see it is not an equilibrium system. We discuss 
the meaning of this observation. 

Finally, we turn to the so-called information paradox. We argue that the fact that the 
geometric optics solution is exact for the two-dimensional black hole tells us that we really 
have to take into account the spacelike line, r = 0, stretching from the quantization surface 
Aq to the surface on which we do measurements. To do this in a consistent way for two and 
four dimensions we introduce a lattice in Lemaitre coordinate rj and demonstrate that the 
spectrum of the time-dependent Hamiltonian is constantly changing; this proves that the 
time dependence of the problem is not a coordinate artifact. We conclude with a discussion of 
the picture which is suggested by this analysis; namely, that a very weakly coupled remnant 
forms as the black hole evaporates. 

VI. THE MOVING MIRROR 

Consider a free field theory in flat space, together with the boundary condition that the 
field vanishes on and to the left of a curve x{t). To discuss the solution of the Heisenberg 
field equations for this problem it is helpful to review the simplest way of solving the free 
field Heisenberg equations when there are no boundary conditions. Start from the free field 
Euler-Lagrange equation 

{d^-dl)<P{t,x)^0, (23) 
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FIG. 3: A plot of the moving mirror and three null geodesies which represent contributions to 
go{x — t). These curves are selected to show that null lines which hit the mirror at later times come 
from points which start out closer and closer to the point x = 1. 



and rewrite it as 



{dt-d,){dt + d,)(P{t,x) = 0. 



(24) 



Then observe that the general solution to this equation can be written as: 



(j)(t, x) = f{x — t) + g{x + t). 



(25) 



The functions / and g are determined by the values of 0(t, x) and its time derivative at 



where 0o and ttq are initial conditions at t = surface. It is a simple matter to integrate 
these equations to determine /(x) and g{x). 



—t + A{l — e"^''*), where we have chosen to plot the curve for A = 1 and k = 1/2. Clearly 
this system has a time- dependent Hamiltonian and, nevertheless, possesses a unitary time 
development operator. If, as shown, we assume that before t = the mirror is at rest and 
has been that way for an infinite amount of time, then it is reasonable to assume that the 
initial state of the problem at t = is the vacuum state for the free field theory defined by 
the condition 0o(O) = 0. 

In this case the Euler- Lagrange equations remain unchanged, however the solution needs 
to be modified to maintain the boundary condition which says that 0(t,x(t)) = 0. This is 



t = 0: 




(26) 



Next consider, as shown in Fig.^, the case of a field theory with moving boundary x{t) 
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easily done by adding a reflected wave go{x — t) to the general solution , so that it becomes: 

(f){t, x) = e{x - t) f{x -t)+ git + x)+ 9{t - x) go{x - t). (27) 

As in the case with no boundary conditions we determine f{x) and g{x) from the initial 
data on the t = surface: 

= ! r<icf^+.(c)l (28) 



2 7o V 

/w-^/:^c(f^-.(o). 

Given this, we flnd from the requirement that the total fleld 0(t, x) vanishes on the curve 
x{t) that go is given by 

go{x -t) = -g{xo{t,x)), (29) 

where, as indicated in Fig.^, xo(t, x) is the point on the t = surface from which the reflected 
wave came. The function Xo(t, x) is determined by observing that the light ray which comes 
to the point (t,x) left the mirror at some point (ti,x(ti)) and that the incident wave which 
arrived at this point came from xo(t,a;). These statements are equivalent to the equations 

X - t = xiti) - ti, xiti) +ti = XQ{t,x), x(ti) = -ti + A(l-e~^'"*i). (30) 

From Eq.(^) we immediately obtain 

a;o(t,x) = A(l-e-2'^*i) (31) 

or, equivalent ly, 

-2t, = ilog(^^^). (32) 

Combining this with Eq.(pOD we derive 

x-t = -2ti + A(l -e"^''*^) = -2ti +xo(t,x), 

x-t = xo{t, x) + i log , (33) 

which, for t — x ^ A, has the approximate solution, 

a;o(t,x) ~ A(l -e-"(*-^)). (34) 

This completes our solution of the Heisenberg equations of motion for the free fleld in 
the presence of a moving mirror. In the sections to follow we will see that because the 
point XQ{t,x) becomes arbitrarily close to the point A for t ^ x, a thermometer placed 
at a distance from the mirror will measure a temperature K/2n and a detector will see an 
outgoing energy flux k.'^/AStt. 
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VII. UNRUH THERMOMETER 



We begin with a precise definition of what we mean by a thermometer. In what follows 
we will take a thermometer to be a simple quantum system, with multiple energy levels, 
locally interacting with the field (j){t,x). Two terms should be added to the massless scalar 
field Lagrangian in order to specify the interaction of the thermometer with the field. The 
first is a term which defines the eigenstates of the non-interacting thermometer and the 
second, is an interaction term of the form 

M„t(t) = ee-(*-*°)'/2'^Q0(a;,t). (35) 

Here e is the small parameter in which we will perturb, to and a define the range in t for 
which the interaction is turned on and x specifies the spatial location of the thermometer. 
The operator Q is assumed to be an operator which causes transitions among the energy 
eigenstates of the thermometer. 

A number of assumptions have to be made in order to get reasonable results. First, in 
order for the thermometer to know the mirror is moving, it is necessary to assume that 
to ^ X- Second, we must impose an adiabatic condition, ^/a' ^ 1/ where E is the typical 
excitation energy of the thermometer, so that we do not excite the thermometer just by 
turning it on or off. Finally, we must impose the condition ~ k, so that the acceleration 
of the mirror is capable of exciting the higher states of the thermometer. 

Given these assumptions, second order perturbation theory in e tells us that the proba- 
bility of the thermometer being excited to a state with energy E is 

r{E,Eo) = e^\{E\Q\Eo)\^ J dtdt'e-'^''-^'>^'''-''^-^^'-''^'+^^^ (36) 

where (O) stays for the vacuum expectation value of the operator O. 

We must now evaluate {(f)(t,x) (j)(t',x)) using Eqs.(p^-|29|), which tell us how to rewrite 



(37) 



(f){t,x) and (f){t',x) in terms of 0(x) and 7c{x) on the surface t = 0. Once we have done this 
we can rewrite these t = operators in terms of annihilation and creation operators and 
evaluate the resulting expression. We obtain 

1 

o2tt (E-Eo)/k _ - 

This result shows that the thermometer reacts as if it is in interaction with a heat bath at 
a temperature k/2tt. From this point on we will assume, without loss of generality, that 
i^o = in order to simplify the equations. 
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^ ' °^ 2 E-Eo 



Let us now discuss the details of the calculation. To obtain Eq. (^) we must first compute 
the quantity 

G{t,t') = {(j){t,x)(l){t',x)), (38) 

when both t and t' are much greater than x. When these inequalities are satisfied we see 
from Eq. ([27|) that (f){t, x) is given by 

<j){t,x)=g{t + x)-g{xo{t,x)). (39) 

It follows from Eq.(^) that 

Using the expansion of the operators 0(C) and 7r((^) defined in Eq. (|2^) we obtain 

dm 

dc 

from which it follows that 



+ 7r(C) = j^du^ (e^-^ ai + e'^^^ a^) (41) 



j roo 1 r / \ / \ n 

(f){t,x) = - duo (e^'^^»(*'^) - e^'^(*+^M al - (e-*^^o(*'^) - e-^^(*+^M . (42) 

2 Jo vvrcij LV /V /J 

It is straightforward to evaluate the expectation value 

{(f)(t,x)(f)(t',x)) = — / — fQi-'^(=^o(*'^)-^o(t,x)) _|_ ^Mt'-t) _ ^iuj{xo{t' ,x)-t-x) _ ^iuj{t'+x-xo{t,x))\ 

' ' 47r i ^ ^ 





(43) 



Computing these integrals we arrive at the result 

{4>{t,x)4>(t' ,x)) = — [\n{xo{t' ,x) — Xo(t,x) + if]) +\n(t' — t + irj) 
An 

— ln(xo(t', x) — t — X + irj) — ln(t' + x — Xo(t, x) + ir])] , (44) 

where an infinitesimal positive quantity r] has been introduced. Substituting this into 
Eq . (^6l) , we see that the last two terms in Eq.(^4|) are damped by a factor of e~'^^^ , which is 
negligible thanks to the adiabatic assumption. The reason this happens is that both t and 
t' are restricted to be near which is assumed to be large enough so that XQ{to,x) ~ A. 
Therefore, either t or t' integration can be done, yielding the suppression factor. This leaves 
only the integration of the first two terms. Since the first term takes the most work let us 
begin with it. 
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The integral we must evaluate to obtain the first term's contribution to the transition 
probability is 

^ r dtdt' e-^^(*-*')e-(*-*o)V2^e-(*'-<o)V2<x x) - xo(t, x) + iii). (45) 

In order to simplify the Gaussian terms it will be convenient to let t ^ t + to and t' ^t' + to 
and then, replacing xq by Eq.(|3TD we obtain 

^ r dtdt' e-^^(*-*')e-*V2<Xg-t'V2<x Z^^^^) _ ^^^^ _ ^ Y^^^-Kt _ g-KtM ^ ^4g) 
47r J-oo ^ ' 

As in the previous discussion, we see that the ln(y4) term and the term — ^(to — x) are both 

suppressed by e^'^^^ . Having removed all of the to dependence it is convenient to define 

u = t' ~t and V = t' + t and rewrite what is left as 

^oo 



, dvdu e^Eu -u^/ia -v^/Aa (v + u) + Xxiie^'' - V 

Stt 7-00 V 2 ^ ^ ^ 



(47) 



Once again all the terms linear in v and u vanish or give a contribution of order e so 
the only integral we have to compute reads: 

(48) 

In order to handle this term we rewrite the integral in Eq. (|i8D as a sum of integrals over 
positive u, to obtain 
-1 



— / du e^^^e"" ln(e"" - 1). 



47r 7o 



e'-^" (ku + ln(l - e"''")) + e"^^" ln(l - e"'^") + 



(49) 



At this juncture we see that if we rewrite the integral in terms of the variable ^ = ku and 
expand the logarithmic terms in a series in e~^, we obtain 

1 / — °° r 

^ n=l 



An K 



oo ^{iE/n—n)^ CO ^—(iE/K+n)^ 



n=l 



n 



(50) 



The condition an? ^ 1 allows us to replace the term e 5^/^'^''^ by unity, which then allows 
us to carry out all of the integrations and obtain 



Ann 



n=l 



2 oo 



+ (f)- 



AE 



(51) 



The sum can be explicitly done using the identity 



E 



n=l W 



1 



+ (f) 



Tin 



^2ttE/k 



K 



(52) 
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which makes the total undamped contribution of this term to the transition matrix element 

^ ^ (53) 



2E 



q2ttE/k _ ^ 

Finally it only remains to show that the contribution of the second term in Eq.(^) vanishes. 
This is easily done by rewriting C, = Eu, using the fact that aE"^ ^ 1 and rewriting 
everything as a sum of integrals over the range u = . . . oo. 

VIII. THE ENERGY FLUX 

Having seen that a thermometer at a fixed location will, at large times in the future, 
measure a temperature T = k,/27t, we would like to delve further into this phenomenon. 
One way to do this is to compute the net flux of energy passing through the point x, to see 
if the thermometer is heating up due to a flux emanating from the location of the mirror. 
This means we need to compute Tf^,, the energy-momentum tensor for the massless field. 

In Minkowski space the energy-momentum tensor for the massless scalar field is defined 
to be 

T^, = ^{d^<p{t,x),d,<P{t,x)}, (54) 

where {A, B} denotes the anti-commutator. In general, the expectation value of any com- 
ponent of Tf^i, is divergent since we have to evaluate the product of two quantum fields at the 
same spacetime point. Thus we need a regularization procedure which can identify possible 
infinities. In what follows we adopt a point splitting method to regularize the stress-energy 
tensor and define: 

Tt. = ^{n{t + 6,x), d^t -6,x)}. (55) 

It is understood that the limit 6^0 should be taken at the end of the calculation. It turns 
out that in this limit the energy density Tu diverges as 1/5^ but the flux, T*^', is finite and 
unique: 

T*^ = — . (56) 

487r ^ ^ 



To derive Eq . (|56|) , one takes appropriate derivatives in Eq. (|42|) and evaluates the expec- 
tation value of the commutator defined in Eq.(|55D The result is: 



T - ^ 



52 (xoiO - xo{^ + 6)y {t + 6 + x-xoiOy {t + x-xo{^ + 6)) 



,(57) 
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where Xg(^) = dxo{^)/d^ and ^ = t — x. Each of the first two terms in Eq. (|57|) diverges if 
6^0 whereas the last two terms take finite hmits and cancel one another. Substituting 
the explicit form of xo{t, x) from Eq. (|3^ into the second term of Eq.(^7|) and expanding it 
in powers of 6 we find that it becomes 

+ (58) 

From this we see that in the limit 5 — > all singular terms in 6 cancel, yielding the result 

T*" = — . (59) 

487r ^ ^ 

There are several features of this calculation which are generic and are therefore worth 
further discussion. 

The first point which merits discussion is the finiteness of our result for the fiux. The 
general reason for that is the theorem discussed in Ref. ^ where divergent terms which can 
appear in T^^, for a general background gravitational field are explicitly given. Evaluation 
of these terms for the case of fiat Minkowski space with a moving boundary as well as the 
Schwarzschild black hole shows that no divergent terms can arise in the computation of 
the fiux and therefore it must come out finite. The preceding discussion shows, by explicit 
calculation, how this works for the case of the moving mirror. 

Next, we wish to observe that Eq.(p7D shows that only those terms which are singular in 
the 5 — > limit contribute to the result for the fiux. A peculiar feature of those terms is 
that for them the points {t + 6, x) and {t, x) can be traced back to the same point on the 
initial surface. Terms which are finite in the limit 6 —>■ always cancel exactly. One way of 
describing this state of affairs is that the first term in Eq. (0) represents fiux coming from the 
right and the second term represents flux coming from the left (i.e., the mirror). In the case 
of Minkowski space with no mirror, these two contributions would cancel exactly, however 
in the case of the moving mirror the flux reflected from the mirror is subjected to a time 
dependent redshift. It is this time dependent redshift which produces a non-cancelling flnite 
addition to the flux. As we will see in a later Section, a similar decomposition of the terms 
contributing to the energy momentum tensor is possible for the case of a black hole. In the 
Schwarzschild case the flux coming from the right behaves much like the Minkowski space 
contribution, but the flux coming from the left emanates from the vicinity of the horizon. 
This flux also sees a time dependent redshift due to the fact that flux which arrives at a 
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slightly later time originates from a point which is slightly closer to the horizon. Just as in 
the moving mirror case, this time dependent redshift is the cause of the finite, non-cancelling 
contribution to the outgoing flux which we identify as Hawking radiation. Finally, we should 
point out that the same approach can be used to describe what happens if the mirror moves 
along an arbitrary trajectory x{t) which asymptotes to the light cone, since all that changes 
in the calculation is the way in which one deals with the function xo{t,x). 

IX. THE BOGOLIUBOV TRANSFORMATION 

Bogoliubov transformations are a common tool used to deal with field theory in curved 
space time, in particular in cases of the moving mirror and the Schwarzschild black hole. In 
this Section we show how to understand these ideas within the Hamiltonian formalism. 

The usual context within which one discusses Bogoliubov transformations is a situation 
in which there is a time ti before which the Hamiltonian is time independent and free, 
and a time t2 after which it is time independent and free. In this case the Bogoliubov 
transformation is straightforward to both define and compute. All one has to do is Fourier 
expand the fields x) and 7r(ti, x) in terms of annihilation and creation operators auj{ti) 
and a^(ti), as in Eq.(^T]) and do the same for the fields (f){t2,x) and TT{t2,x). Then, if we 
have explicit formulas relating the fields (f){t2,x) and 7r(t2,a;) to the fields at time ti we 
get an explicit relation between and a^^lti) and their counterparts at time t2- This 

relationship is the desired Bogoliubov transformation. Applying these ideas to the case of 
the moving mirror would be quite simple except for the fact that, although before time 
ti = the theory is indeed time independent and free, there is no time t2 for which the 
same conditions apply. There is, however, a sense in which we can treat the system as being 
almost time independent. We can then use the corresponding Bogoliubov transformation to 
compute such things as the thermometer response and outgoing energy fiux. 

To understand why this works consider the coordinate system shown in Fig.^ and focus 
on the point {t',x'). Causality requires that for times t such that x + 1' > t > t' (i.e., for 
all points {t, x) lying within the shaded triangular wedge) the field operators x') and 
TT{t',x') evolve according to the infinite volume free field equations. Thus, to compute what 
a thermometer moving along the timeline {t, x') would see, we can expand the initial t = 
state in terms of the annihilation and creation operators of the instantaneous Hamiltonian 
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FIG. 4: The coordinate system used to discuss the BogoUubov transformation. Note that before 
time t = and in the shaded wedge formed by the hues t = t' and the 45° hne we use ordinary 
rectangular Cartesian coordinates, but in the remaining region we use hues which are translations 
of the mirror trajectory. For fixed x the field (/>(r, x) for all (r, x) in the shaded wedge can be 
computed in terms of its values on the line t = t' using free-field equations of motion. 

H{t') and then use infinite volume free field equations to evaluate the relevant correlation 
functions. This approach should lead to our earlier results if, as we will show in a moment, 
we work at large enough x. It has (obvious) limitations, however, if we want to study this 
problem at some fixed finite distance from the origin over a finite interval in time. 

From Fig.H we see that the farther out x' is along the x-axis the longer the time interval 
for which infinite volume free field propagation makes sense; nevertheless, we see that for any 
finite x' the approximation eventually breaks down. This means that Fourier transforming 
in the time variable in order to identify the annihilation and creation operators is only an 
approximation to what is going on. One always has the option of putting the initial surface 
in the infinite past and working at x' = oo, which would make the approximation exact; but 
if one is interested in a detailed picture of how the system develops in time, working at past 
and future infinity is a severe limitation. Having said that, let us show how the Bogoliubov 
transformation should be understood in the Hamiltonian framework. 

To carry out the Bogoliubov transformation we begin by considering the instantaneous 
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Hamiltonian of the system at time t = t', when the mirror is at the position x{t'): 

H{t') = - dx (nit', xf + d^cPit', xY 
2 Jt' ^ 



(60) 



This is just the Hamiltonian of a free massless field theory defined on the interval x{t') < 
X < oo and it is diagonalized by expanding the field and its conjugate momentum in terms 
of annihilation and creation operators and bl., defined by the formulas 



°° dk 



V/evt 



sm{k{x - x{t'))) (bl + bk 



Tr{t',x) = —i J dk]j— sm{k{x — x{t'))) (bl — bi, 



(61) 



(62) 



It will be convenient in what follows to make the coefficients of bl and b^ simple expo- 
nentials by taking the linear combinations 



TT{t',x) 



-7r{t',x) + 



d<p{t', x) 

dx 
d(f){t', x) 



dkJ- (bl e*'^(---(*')) + bk e-*'^(---(*'))) . 



(63) 
(64) 



dx JO V TT 

The next step in computing the Bogoliubov transformation is to rewrite </)(t', x) and 
7r(t',x) in terms of the t = fields. Following our earlier discussion we see that 4>(t,x) is a 
sum of two terms: 



0(t', x) = e{t' - x) (j)i{t', x) + e{x - t') (j)ii{t', x). 
and, using Eq.(H^), we obtain the following expressions for 



7rjit',x) + -^{t',x) 



dk 



Ik r 



and 



,,(t',a;)-^(t',x) = ^^°^''"^ 



dk \ l — 

TT 



aie 



ikxo{t'x) 



+ ak e 



-ikxo{t' x) 



(65) 



(66) 



(67) 



dx dx 

If we consider the last of these two equations and ignore the part of 0/(t', x) which does 
not come from the mirror, and if we assume that t' is large so that x{t') is large and negative, 
then to a good approximation we can extract the coefficient of bl by simply taking the Fourier 
transform 

' , + r dx 



■ik{x-x(t')) 



1Tl{t',x) - — {t',x) 



k(x-x(t')) ^^o(J^''^) ^^^iiJxo{t',x) _^ ^^^-iojxoit' ,x) 
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Changing variables from x to xo{t', x), rewriting Eq.(|33|) in the hmit of large t' as 

x = t' + - ln(^-^) + xo (69) 
and letting xq = ^(1 — ^^e expression for b\ becomes 



71 J Zn \ 71 Jo 

Since the dominant contributions to this integral come from small ^ we can extend the 
^-integration to infinity and evaluate the result in terms of Gamma functions to obtain: 



(71) 



A{uj-k)J \A{uj + k)^ 

Given this approximate Bogoliubov transformation one usually computes the expectation 
values of the number operators (0| 6^,6^10), where (0| stands for the state the system started 
in at time t = 0, which we have chosen to be the state annihilated by the operators a^. 
Using Eq.(|7l| ), one derives 

uduj 



(72) 



Except for the divergent integral in front of the expression, this is what we need to plug into 
the formula for the transition probability for the thermometer to obtain its response. The 
usual way of arguing away the divergent prefactor is that the number of particles grows with 
time and that one should divide by this logarithmically divergent term to get a number of 
particles per unit time. If one does this then one can use the results of this approximate 
calculation to derive the response of a thermometer which is switched on and off for a finite 
period, or the flux of energy through a given point. 

To conclude this Section we would like to point out that we do not see any advantage, 
either computational or conceptual, in discussing particular choices of the vacuum states 
and the associated Bogoliubov transformations for this explicitly time dependent problem. 
The time dependent nature of the problem, which explicitly manifests itself in the time 
dependence of the Hamiltonian, naturally leads to a long time steady state behavior of 
physical system that is practically independent of the initial conditions. 
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X. TWO MIRRORS: A CURIOUS PHENOMENON 



Let us complete our discussion of flat-space problems by considering a variant of the 
moving mirror in order to show that for this problem, at late times, almost all of space- 
time will be filled by fields that correspond to degrees of freedom of the massless theory 
which were originally concentrated in an exponentially small region of a single point. This 
peculiar, but otherwise absolutely trivial property of the massless theory in the presence of 
a time dependent Hamiltonian, will be encountered again when we discuss the question of 
Hawking-Bekenstein entropy for the case of the black hole. The only modification we will 
make of the original moving mirror problem will be to add a stationary mirror at x = L, see 
Fig.|. 

A glance at Fig.^ shows that adding the second mirror does not complicate the compu- 
tation of the field operators at any time in the future, since at most two reflections should 
be taken into account. The interesting feature of this solution to the field equations is that 
now, for times t > L and points (t, x) lying to the right of the shaded region bounded by 
the mirror and the line x = 2L — t, there are no direct contributions to the field operators. 
Instead, the value of the fields at such a point is the sum of contributions coming either from 
a single reflection from the moving mirror or a reflection from the moving mirror followed 
by a reflection from the stationary mirror. Thus, we see that for large t and points {t, x) 
lying to the right of the shaded region, the fields 4>(t,x) and n(t,x) only depend upon the 
fields on the original t = surface which lie within an exponentially small neighborhood of 
the point x = A. 

XI. THE BLACK HOLE: GEOMETRIC OPTICS APPROXIMATION 

Finally, let us turn to a discussion of the Schwarzschild black hole. We pointed out in 
Section IV that by introducing Lemaitre coordinates we can canonically quantize the theory 
of a massless scalar field and use the resulting time-dependent Hamiltonian to derive the 
corresponding Heisenberg equations of motion, Eq. (|T9|) . As in the case of the moving mirror, 
the behavior of the system at later times is obtained by solving for the Heisenberg fields 
as a function of the fields defined on the initial surface on which we quantized the system. 
In the flat space case we argued that the easiest way to solve the Heisenberg equation of 
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FIG. 5: A plot of null geodesies which contribute to go{x — t) at late times in a problem with one 
moving mirror and one fixed mirror. Note that eventually the only contributions to fields which 
lie outside the shaded region come from geodesies which reflect once from each mirror or geodesies 
which reflect only from the moving mirror. 

motion is to trace back the two null-rays leaving the point {t, x) (the point at which we wish 
to compute the field and its conjugate momentum), to the two points from which they leave 
the t = surface and write the answer in terms of the and vr at those two points. We will 
now show how to generalize this approach to the corresponding equations in curved space. 

To motivate what we call the geometric optics approximation let us study Eq. ([T9D in 
Painleve coordinates (A,r), since these coordinates are non-singular and the dependence of 
the solutions on A and r factorizes. The usual WKB approach to this sort of problem is to 
assume a solution of the form 0o = T~^e*'^'^/t^(r) and substitute this ansatz into the field 
equation. In this way one obtains that, for large uj, fuj{r) can be written as 



We now observe that these solutions are constant along incoming or outgoing null geodesies 
where an incoming null-geodesic starting at the point Xi at time A = is a curve r(A) such 
that 



and similarly, an outgoing geodesic starting at X2 at A = 0, is a curve r(A) such that 



ln/c^(r) 



^uJSl,2{r) + 0{u-^), Si,2(r) = ±r - 2v^ ± ln((v^ ± l)^). 



(73) 




(74) 



S2{x2) = X + S2{r{X)), 



(75) 



where 5*1^2 are as defined in Eq. ([73D . 
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To change this observation into an ansatz for the solution to the S-wave field equation, 
we simply mimic the general form of the solution for the moving mirror in fiat space; i.e. we 
say that for general (A, r) 

0o(A,r) = l(0i(A + ^i(r)) + <^2(A + ^2(r))), (76) 

and the functions (t>i,2{Si^2{r)) = fi,2{r) are to be determined from the boundary conditions 

0o(O,r) = ^, a,0(A,r)U=o = v^^i(r), (77) 
r 

where 0i(r) and 7ri(r) are the rescaled operators we introduced to quantize the theory on 
the initial surface A = 0. 

Substituting Eq.([75|) into Eq.(^) and following the same procedure as in the fiat space 
case we obtain 



(78) 



^3/2 
u 

where (f)[ = d(f)i/dC, and 5*1^2(2^1,2) = A + Si^2{r). Combined with the fact that field (pi and 
its momentum tti are expressed through the creation and annihilation operators defined at 
A = 0, the above set of equations allows us to compute any Green's function of the field 0o 
at any later time. 

Let us note that for the four-dimensional theory the geometric optics form of the solution 
is not valid as r 0. This is a point we will return to when we discuss the information 
paradox in the last section of this paper. 



XII. BLACK HOLE: THERMOMETER REDUX 



As for the moving mirror, we expect the time dependence of the Hamiltonian to refiect it- 
self in the existence of steady state phenomena such as outgoing radiation. We will now show 
that this is indeed what happens and that this formalism naturally leads to the prediction 
that there is an outgoing, time-independent fiux of energy which at large distances corre- 
sponds to a body at a Hawking temperature Th. To do this we follow the procedure used for 
the moving mirror and weakly couple the massless field to a detector [Q, |^ (which acts as a 
thermometer) located at some fixed Schwarzschild (or Painleve) radius r. As before, we add 
an interaction term to the free field Lagrangian of the form M^t ~ e-^^-^<^'^'/^'^^^'^Mt,r)Q, 
where Q is an operator which acts in the Hilbert space of detector eigenstates. Second order 
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perturbation theory in V^nt tells us that the probability of exciting the detector to a state of 
energy E is related to the Fourier transform of the Green's function of the massless field [0: 

V{AE) ~ \{E\Q\Eo)\^ J dMt'e-'^^(*-*')-[(*-*°)'+(*'-*°)'l/(2''')(0o(t,r)0o(t',r)), (79) 

where AE = E — Eq and Eq is the ground state energy of the detector. The only difference 
is that now the Green's function in Eq. (|79D is to be computed using the evolution equation 
for the field (j)o{X,r) which relates it to 0i(r) and 7ri(r) on the surface A = 0. 

As in the case of the moving mirror it is convenient to define the points xi^2 and yi^2 as 
follows: 

Si{xi) = Ai + S'i(ri), 5'2(x2) = Ai + S'2(ri), 

Siivi) = A2 + S^{r2), = X2 + 52(r2), (80) 

so that Xi and X2 are the points on the A = surface from which infalling and outgoing 
null geodesies must leave to arrive at the point (Ai,ri) and yi and ?/2 are the points from 
which infalling and outgoing null geodesies must leave the same surface to arrive at the 
point (A2,'r2). Since we assume that the thermometer stays fixed at the same Schwarzschild 
r we must identify ri = r2 = r and, also, we must remember that there is a transformation 
between A and t so that Ai = Ai (t,r) and A2 = A2(t',T). At this point we need to evaluate 
Xi^2 and yi^2 for large values of to- Given the explicit form of the functions S'1^2 it follows 
directly that these equations have the approximate solutions: 

= t + r, X2 = 1 + 2e~(*^'')/2^ 
2/1 = t' + r, ^2 = 1 + 2e-(*'-^)/^ (81) 

where we have assumed that t t' to- Asymptotically, X2 y2 —>■ ^ and xi ^ yi —>■ cxd. 
In this limit, the Green's function can be written as 

— 1 / in \ 

{(j)o{t,r)(f)o{t',r)) ^ ^— ^ Mn \xi - yi\ + ln\x2 - 2/2I + y [K{xi,yi) + K{y2,X2)] + cj{^2) 

where k{x, y) = 6{x — y) — 0{y — x) and c is some constant. 

It is instructive to consider the terms in Eq.(^) separately. The constant does not 
contribute to V{E) since, as in the case of the moving mirror, it yields a result proportional 
to exp{—AE'^6^) <^ 1. The In \xi — yi\ term and the terms described by the function k give 
simple contributions that can be written as 

V,{AE)^-^ + 0{i5AE)-'). (83) 
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The important part of the final result comes from the second term in Eq. ([8^ ) which describes 
the radiation coming from the vicinity of the horizon. Appropriately shifting the integration 
variables and restoring factors of 2M we obtain 



If we then change the variables to v = {t + f) , u = (t — t') and neglect all the suppressed 
terms we arrive at 

OttS r 1 n 

(85) 



1 l' 

gAS/TH -1 + 2. 



where the Hawking temperature Th = l/(87rM) has been introduced. Since the total prob- 
ability is given by the sum of Vi and V2 we have the final result: 

n^E) mQ\E,)\^ 1 .... 

—5 ^ ^KeIt^i (86) 

The interpretation of this formula is straightforward. If, at a large, fixed distance from the 
black hole, an observer switches on a detector which interacts with the massless field for 
finite amount of time, then the energy levels of the detector get populated as if the detector 
was in equilibrium with a thermal distribution of particles at a temperature Th = l/87rM. 

The only subtlety which we have skipped over in this calculation done for large r, is that 
for arbitrary r the interaction term should have a correction for the time dilation at point r, 
since the energy levels of the thermometer are defined in its rest frame. The result of adding 
this to the calculation is that, putting back factors of M, a thermometer at arbitrary r will 
read a temperature 

kBT = , ^ (87) 

MM^l - 2M/r 

XIII. BLACK HOLE: ENERGY FLUX 

The Schwarzschild calculation of the energy fiux through a sphere of fixed radius is done in 
the same way as for the moving mirror: we point-split the fields appearing in the expression 
for the energy momentum tensor, regulate the resulting expression and then take the limit 
of zero splitting. The result of this computation is that we find that Tx ^, is finite and 
non- vanishing and the total fiux through a sphere of large radius is given by 

Flux = ^ = (88) 
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The full expression for Tx^r) contains terms which vanish for large A and therefore can be 
identified as transients; persistent terms which decrease faster than and, finally, per- 
sistent terms which fall off as and hence contribute to the flux. While we have carried 
out the computation for arbitrary (A,r), the resulting expressions are too cumbersome to 
present here and we limit our discussion to the computation of terms that both approach a 
constant for large A and fall off as 

Let 0" = A2 — Ai and define the flux as: 

(Ta,) = hm -(0| [ ^ \ |0)- (89) 

The formula for the total energy passing through a large sphere, in Lemaitre coordinates, is 
given by 

r lT^ ) 

J = Ita / d^d»^9"9'"^ = - Hm p^(T„>, (90) 

where a normalization factor of (47r)^^ is introduced since (po denotes the S'-wave component 
of the massless field. Given this expression, it is straightforward to compute this flux using 
the explicit expression for the time evolution of the field (po, Eq-(^), and then take the limit 
A2 — ^ Ai. 

Since in the geometric optics approximation the field and its conjugate momentum are 
given in purely geometrical terms (i.e., in terms of the point at which the field is to be 
evaluated, the null geodesies arriving at that point and the two points on the initial surface 
of quantization from which they left), we are free to carry out the computation in any 
coordinate system. Since the calculation of the flux for large values of A and 1] is simplest in 
Painleve coordinates we will transform to these coordinates in the discussion which follows. 
However, because we define the point-split energy momentum tensor by holding rj fixed and 
separating the fields in the time variable A we have to take into account that the variables 
ri and r2 corresponding to (Ai,//) and {X2,v) ^-re related by 



= - ^ - ^ - -^-i^ + (^1) 

Taking the derivatives in Eq.(|89D and considering the limit ri 2 00, we obtain the following 
expression for the flux 

J = -l lim [({/{ (xi), f[{yi)})W,{r2, ri, x,,y^) + ({/^(xa), /^(y2)})W^2(r2, r,, x^, I/2)] , 

(92) 
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where the functions /i,2 are defined in Eg. (ffSD , and 



W,{r,, n, X, y) = ^TTY^ fl - ^1 • (93) 

S'i[x)S[[y) \ J¥{ 



Once again (cf. Eq.(pOD) we define xi,2, 1/1,2 to be the points on the A = surface from 
which the null-geodesics which wind up at the point (A, 77) originate. Given these equations 
it is easy to derive the relation between 7/1^2 and xi^2 as a power series expansion in a. It 
follows from the specific form of the functions 5*1^2 that, in the limit of large A and large rj, 
the limiting values for these points are xi — 00, X2/2M 1. 

To compute the expectation value of the anticommutator ({/{(xi), /((yi)}) we start from 
Eq.([78D and derive 

{{f[M. f[{yi)}) = \ [({vri(xi), My,)}) + <P[{yi)}) + ...], 

where the dots represent terms which do not contribute to the fiux in the ri — 00 limit and 
the anticommutators are computed using Eq.(^). In this way we find 

m{x^),<P'M}) = v?^^^' -1(1/1)}) = -;77^^' (94) 

so that the final result for the anticommutator reads 

({/{(^i), f[iyi)}) = ^ ,2 + - (95) 

2n [xi - yi)^ 



Performing a similar calculation for the second term in Eq. (|92D and substituting expansions 
of yi^2 and r2 in terms of 0:1,2 and ri, we find for the energy fiux 

1 TT 

^ " 1927r(2M)2 = 12^^' ^^^^ 

where once again we have introduced the Hawking temperature Th = l/(87rM). 

Eq.(PBD shows that the energy fiux at large distances is finite. We have already noted 
that this result is in accord with a general theorem that deals with the structure of the 
possible divergences in the stress-energy tensor computed in a gravitational background . 
It turns out that one can derive an explicit formula for all of the possible divergences which 
can occur as coefficients of specific functions of the metric and renormalize them by adding 
explicit counterterms to the Einstein Lagrangian. Adding these terms to the Lagrangian 
and computing the resulting modifications of the energy momentum tensor, one finds that 
there are no terms which can contribute to the ojf- diagonal element of the energy momentum 
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tensor in both the flat space and the Schwarzschild metric in Lemaitre coordinates. Thus, 
since there are no possible counterterms which can remove divergences in the flux, the result 
must come out finite as it indeed does. Unfortunately, for finite values of r our result still 
contains logarithmically divergent terms 0{lna) multiplied by functions of r that decrease 
faster than l/r^ in the limit of large r. We believe this to be due to the fact that we 
have restricted attention to the L = component of the field and have not considered 
higher angular momenta. This observation is supported by the fact that in the case of a 
two-dimensional black hole, where higher angular momentum modes are absent, our result 
for the flux is finite for all values of r. 

To conclude our discussion of the flux let us comment a bit on the back reaction issue. The 
problem of back reaction is equivalent to the statement that the computation we are doing 
is not, from the point of view of the Einstein equations, self-consistent even at the semi- 
classical level since, on the one hand, for a static Schwarzschild metric the Einstein tensor 
Gf^i, vanishes for all r 7^ but, on the other hand, we find that the energy-momentum tensor 
of the scalar field has, at large times, a finite, uniquely defined, off-diagonal component in 
this background. 

Since our approach studies the behavior of the scalar field theory starting from a well 
defined quantum state at a finite initial time, we should be able to solve the problem itera- 
tively, computing corrections to the Schwarzschild metric due to non-zero expectation value 
of the stress-energy tensor at the right hand side of the Einstein equations. Although it is 
quite difficult to do this in general, we would like to point out that a simple modification of 
the metric: 

ds^ = _(l - ^M^) dt' + \^ dr' + r'd^\ (97) 

produces the Einstein tensor of the form 
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(98) 



where X(t) = 1 — 2M(t)/r and where we only retained terms which are linear in dM(t)/dt. 

Although the approximation for the metric Eq. ( pT]) is too crude to reproduce the vacuum 
expectation value of the stress energy tensor of the massless field for arbitrary values of t and 
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r, for large values of r we can easily match the expression for the Einstein tensor, Eg. (pSD , to 
the vacuum expectation value for the energy flux. We then obtain the well known equation 
that describes the evaporation of the black hole: 

^ = M V (99) 
dt 12 Wm) ■ ^ ^ 

In principle, it should be possible to iteratively improve on this approximation and in 
this case solve the back reaction problem at the quasi-classical level until relatively late in 
evaporation process, thus obtaining a better insight into the black hole dynamics. 



XIV. AN INFALLING REFLECTING MIRROR 



We have now discussed the canonical quantization of a massless scalar field in the back- 
ground of an eternal Schwarzschild black hole and have shown that, in the geometric optics 
approximation, one can compute all of the usual results such as the Hawking temperature, 
as measured by a thermometer held at fixed Schwarzschild r, and the flux through a sphere 
of the same radius. Now we would like to address the question of the choice of initial state 
on the surface of quantization. We have already demonstrated that on any fixed Painleve or 
Lemaitre time-slice the Hamiltonian exhibits no singularities and, in fact, is the Hamiltonian 
of an ordinary massless free field theory. Thus, there is no problem in defining the A = Aq 
Hilbert space. 

The only issue which remains is a particular choice of the initial state for the field theory 
on the quantization surface. At this point we would like to stress that practically for any 
choice of the initial state at A = surface which differs from the vacuum state by a finite 
amount of energy, the large time behavior of the system will be exactly the same as described 
in the previous Sections. However, we should also note that a sufficiently bizarre choice of 
the initial state can result in no Hawking radiation at late times. 

To show that such choices of the initial state are unnatural, we would like to pose a 
problem which uses a stationary reflecting mirror to guarantee that for inflnite times in the 
past there is a natural choice of initial state. Next, we let the mirror start to fall into the 
Schwarzschild black hole along a Lemaitre timeline and ask if, starting from this well deflned 
initial state, we see Hawking radiation at late times. We will argue that we do and that this 
easily follows from the formalism that we have already discussed in detail. 
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FIG. 6: The curved dark line represents a mirror which initially, for infinite times in the past, sits 
at Schwarzschild r = 1.4 and then starts to fall into the horizon along a trajectory r(A, r/). The 
horizontal and vertical lines represent an incoming null- geodesic hitting the infalling mirror just 
before it crosses the horizon and being reflected outwards. From the picture it is clear that the 
closer it is to the horizon when it is reflected, the longer it takes to get out to large Schwarschild r. 

The problem we have in mind is a generalization of the problem of the moving mirror to 
the case in which the mirror moves in the background of a Schwarzschild black hole. To be 
specific we assume, as is shown in Fig.^, that we are dealing with a perfectly reflecting spher- 
ical mirror (i.e. a surface on which we assume the scalar field must vanish) of Schwarzschild 
radius 1. Furthermore, as with the moving mirror we assume that the field 0(r) exists 
only in the region outside the mirror. Next, we assume that at some finite Painleve time Ai 
this mirror starts to fall in along one of the Lemaitre curves r(A) = (i?^/^ — 3(A — Ai)/2)^/^. In 
the geometric optics approximation this problem is solved in essentially the same way as in 
our earlier discussions, except now we have to take into account contributions coming from 
null geodesies which fall in towards the mirror and then are reflected out. If we quantize 
on an initial surface Aq <^ Ai, then for all times less than Ai the field in the region r > R 
comes from geodesies which propagate in, essentially, flat space. Comparing this with our 
previous calculation we see that for this region there is no generation of Hawking radiation. 
Remember, that in order for the thermometer to measure a Hawking temperature, or for 
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the energy momentum tensor to show a non- vanishing flux of apparently thermal radiation, 
two events which are separated by a small interval in A must, when traced back to the 
initial surface of quantization by null geodesies, come from points which are exponentially 
closer together than when they started out. We saw in the previous discussion that this 
only happens when these null geodesies pass very close to the horizon. This means that the 
radiation comes from the rays which are reflected from the mirror as it passes through an 
exponentially small region near the horizon (see Fig.|^). From this argument we see that the 
Hawking radiation seen at large future times is generated from points on the initial surface 
of quantization which all lie in a small neighborhood of the point ro, defined by the condition 
that an infalling null geodesic drawn from this point would hit the mirror just as it passed 
through the horizon. This point tq plays essentially the same role as the corresponding point 
Xo{t, x) in our original discussion of the moving mirror. 

Therefore, in order to suppress the Hawking radiation in this setting, one would need to 
significantly modify the vacuum expectation values of field operators and their products at 
what is essentially an arbitrary point, tq. It is this strange requirement on the initial state 
that we consider to be unnatural and from this we conclude that the Hawking radiation is 
a robust phenomenon practically independent of the initial state. 

XV. ENTROPY 

The next topic we would like to discuss is the question of the entropy of the black hole. 
Clearly if we deal with a Hamiltonian system that starts in a pure state and experiences 
completely unitary evolution it has zero entropy (since the entropy of a pure state is zero) 
despite the fact that it exhibits all of the phenomena associated with Hawking radiation. 
Nevertheless, despite this obvious statement we can slightly modify the problem we have 
been discussing and construct an object which, from the point of an outside observer appears 
to have an energy M, temperature = l/SvrM and an entropy 5* = Ah/ 4:, where is 
the area of the horizon. Ah = AirM"^. Let us see how this is done and then ask what is 
happening. 

For this purpose consider a black hole surrounded by a perfectly reflecting mirror of a 
large radius R. In contrast to the infalling mirror case, however, we now assume that the field 
degrees of freedom live inside the spherical mirror; i.e., in the region < r < R. By putting 
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this mirror around the hole we don't allow any energy to escape the region surrounded by 
the mirror, since T^^ is locally conserved; hence the total energy of the enclosed region 
is, for all times, given by the original mass of the black hole M, which can be measured 
by an outside observer, simply by dropping a test particle and measuring its acceleration. 
Next, imagine that the outside observer has a thermometer inserted through a very small 
hole in the mirror. Clearly, nothing is changed in the calculation presented previously and 
so this thermometer measures a Hawking temperature Th- Thus, we see that we have a 
macroscopic object which appears to have an equation of state that says the energy of the 
object is inversely proportional to its temperature. If we now try to use thermodynamic 
concepts, the usual formula dU = T dS becomes 

dM = dS, (100) 

so that by integrating it and by setting the entropy of a zero mass object of this type to 
zero, one derives: 

S = AttM'' = (101) 

This is the usual Bekenstein argument. Note however that the entire discussion wc have just 
given is being applied to a quantum system during the period when it is certainly described 
by a pure state. This makes it problematic to assign a non- vanishing entropy to the system 
and so are we are left asking the question what the entropy of the system constructed in 
such a way means and what it tells us about the system. 

In contrast to the standard notion of entropy in the equilibrium thermodynamics, the 
entropy of the system black hole plus the massless field docs not tell us anything important 
since we must only look inside our reflecting mirror and note that the system whose proper- 
ties are being measured is never in equilibrium. This lack of equilibrium is not because the 
black hole is evaporating, but because the apparently thermal flux arriving at the mirror 
comes from the horizon, is reflected and then disappears through the horizon to be stored 
at r = 0. What we arc looking at is a steady state phenomenon in which an exponentially 
small region near the horizon serves as a constant source of new radiation with essentially 
the same properties as the reflected radiation which then disappears behind the horizon. At 
flrst these assertions might seem peculiar, but we would remind the reader that this sort of 
phenomenon was already seen in the case of the two-mirror problem. Another feature of the 
two- mirror problem which is intimately related to this behavior is that the two null geodesies 
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attaching to the point (A, r), for large enough A, either come directly from an exponentially 
thin region around the horizon, or are rays which originated from such a region at an earlier 
time and then were reflected back to arrive at (A, r). Hence, in the geometric optics approx- 
imation, we see that for sufficiently large A the fields between the horizon and the mirror 
are functions of fields which lay within an exponentially small region around the horizon 
on the initial surface of quantization. This phenomenon is the source of both the Hawking 
radiation and the reason why late time Green's functions seem to show no correlations over 
finite separations in time and space. Prom this argument we see that the Bekenstein entropy 
is somehow related to the curious property seen in the problem of the moving mirror, that 
the causal structure of the theory guarantees that after a finite time the fields between the 
horizon and rcfiecting shell are only functions of the degrees of freedom localized within an 
exponentially thin shell surrounding the Schwarzschild radius. At best, for this problem, 
the Hawking-Bekenstein entropy is a reflection of this fact. 



XVI. THE INFORMATION PARADOX 



Given that a Hamiltonian formulation of the problem of a massless scalar fleld in the 
background of a large Schwarzschild black hole exists, there cannot be an information para- 
dox until one comes to grips with the question of what happens during the final moments as 
the black hole evaporates. While our approach doesn't allow us to discuss these violent final 
moments of the process of black hole evaporation, it does provide insight into the question 
of where and how degrees of freedom which "fall into the black hole" are stored. It also 
provides a different picture as to what might happen after evaporation has taken place. It 
is this picture of the evolution of the problem from its initial state to late times that we 
discuss in this section. 

As noted earlier, in the case of a four-dimensional black hole the "geometric optics" 
approximation fails to give an accurate description of what is happening near r — but, 
fortunately, this failure of the approximation is not an insuperable obstacle to obtaining a 
more complete understanding of the physics. There are two reasons for this. First, if we 
restrict attention to the case of a two-dimensional black hole (i.e., the theory obtained if 
we restrict the metric to just the upper 2x2 matrix g^,^, /x, z/ = 0, 1), the geometric optics 
approximation is in fact exact. Second, there exists a useful discretization of the problem in 
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Lemaitre coordinates which allows one to consistently investigate the problem in both two 
and four dimensions. 

The importance of the fact that the geometric optics approximation is exact for the case of 
the two dimensional black hole is that it tells us what is happening at r = 0. Our treatment 
begins by canonically quantizing the massless scalar field theory on a surface of constant 
Painleve time, A = Aq, and so we are free to require that on this surface the field vanishes 
at r = 0. For subsequent times, however, this boundary condition will be true if and only 
if it is consistent with the Heisenberg equations of motion and a simple computation shows 
that for A > Aq the field does not vanish. 

To understand the context of this observation, let us note that the line (A, 0) is spacelike. 
Thus, what we learn from the solution to the field equations is that when we restricted the 
integration over t] to run from X < rj < oo, in Eg . {\L^ , we were not doing the correct thing; 
we should have included the spacelike surface r = running from Aq to A. This prescription 
should have been obvious from the situation shown in Fig.^, where we see that as A increases 
the surface of fixed Painleve time Aq gets mapped onto the line r = and the surface of 
fixed Painleve time A. Given that we know this prescription is required in two dimensions 
it is not much of a stretch to assume that the same is true in the four dimensional case. 

Because, in four dimensions, the geometric optics approximation is not valid near r = 
we need to do something else to get a better feeling for the problem, something which agrees 
with the geometric optics approximation in two dimensions. To arrive at such a treatment 
of the problem we need to deal with two issues. The first issue has to do with the fact that 
r = is the location of a true singularity in the metric, where the curvature diverges and 
one expects quantum gravity to play a role. Thus a full treatment of this problem would 
perforce need to go beyond the semi-classical approximation. A less profound technical 
problem is that the Hamiltonian treatment along the line r = needs to be carefully done, 
especially at the points where the line r = and the surface of constant Painleve time meet. 
While we have nothing to say about what the true quantum completion of the theory of 
gravity might be, the second problem is easily handled if we assume a minimum value for 
'"min = e, formulate the Hamiltonian problem and then take the limit e ^ 0. For the two- 
dimensional problem this amounts to quantizing the field theory on the surface of Painleve 
time Aq imposing the condition that the field vanishes at r = e and then using the geometric 
optics construction to evaluate the field and its conjugate momentum along the curve (A, e) 
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in the future. Unfortunately, this simple construction is not sufficient to handle the four 
dimensional theory. 

In order to study both the two and four dimensional problems in a consistent manner, 
we propose to discretize the problem by introducing a lattice in rj. To do that we replace 
the continuous variable by a discrete variable rjj defined by 

Vj = S{e + j) (102) 

where 5 is a lattice spacing which has dimensions of {lengtKf/'^ and e is a dimensionless 
parameter such that rmin = (35e/2)^/^. With this definition of the lattice we then introduce 

dimensionless rcscalcd fields $, 11 and the rescaled time variable 

$(A, 77) = 5^/'ct>{\ 77), n(A, 77) = 5'/M\ v), A ^ ^ (103) 



and rewrite the time dependent Hamiltonian as: 

"2n(A,j)' /3\5/3 



.7 = 1 



+ (2) v{\jf'm,j + ^)-H\j)f 



(104) 



3r7(A,j) 

using an additional assumption that 

77(A, j) = ee{X-e~j) + {e + j- A) e{e + j - A). (105) 

This is in accord with our previous remarks which say that once the field reaches the point 
77(A) = 6{e + j — X) = Se it stays there. 

At this point we should make two observations. First, it is obvious that the number of 
degrees of freedom remain the same in this latticized version of the massless field theory and 
the entire effect of the metric appears in the time dependence of the coefficients appearing in 
the Hamiltonian. The second, is that the lattice we have introduced is somewhat peculiar. 
This is because spacing lattice sites equally in the variable 77 does not correspond to spacing 
them equally in Schwarzschild r. In fact, since r(A) = (3(r/ — A)/2)^/^ we see that for 
large values of r the spacing between two neighboring lattice points decreases like S/y/r. 
Thus, while the lattice provides a good cutoff for the field theory inside the Schwarzschild 
radius it is not an effective cutoff for the theory at large r and the discussion we gave for 
continuum theory at large t and r continues to be necessary. Since, for now, we are most 
interested in the behavior of the theory near r — 0, this latticized version of the field theory 
is what we need in order to discuss the physics in a way which goes beyond the geometric 
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optics approximation. Because the lattice Hamiltonian we have introduced is exphcitly time 
dependent, a full discussion of the physics of the theory would require computing the unitary 
time development operator U{t), which is beyond the scope of the present paper. Despite 
this, we can use this formulation of the theory to gain some better understanding of what is 
happening and to discover a peculiar property of the combined system of the massless scalar 
field in a Schwarzschild background. 

The first question which we can address using this formalism is whether the apparent time 
dependence of the problem is a coordinate artifact, in that we have chosen a coordinate sys- 
tem which, although free of singularities, makes the metric time dependent. Fortunately, 
we can buttress our claim that the problem is intrinsically time dependent, even without 
computing the operator U{t), by observing that the spectra of the instantaneous Hamilto- 
nians H{X) are changing as a function of time. Qualitatively we can see that this is the case 
by comparing the spectrum of the Hamiltonian, Eq. ( |104|) , for A = 0, against what it would 
be for a relatively large value of A. When A = and 6 is very small this Hamiltonian will 
be a discrete version of the continuum Hamiltonian we discussed earlier and we expect the 
spectrum to be that of the zero angular momentum mode of a free massless field theory, i.e., 
proportional to k'^. Now, when A is large, then approximately A lattice sites lie on the curve 
r = e and the part of the Hamiltonian H{X) which refers only to these sites has the form 

Hi = Y: + i-ef/' (0(A, J + 1) - 0(A, j))^ (106) 

which after rescaling H(j) to absorb the factors of e, becomes 

= E ^ + ^4 (|j (0'(A,j + 1) - 0'(A,j))^ (107) 

Since this has the form of a latticized nearest neighbor interaction we see that for large A 
the spectrum of the kinetic term will be 

£(A:)=4(-j (l-cos(fc)) (108) 

where < /c < 27r/A since A is the number of sites on the lattice. The remaining part of 
the Hamiltonian consists of two pieces. The first piece, which we will ignore for now, is a 
single term linking the curve r = e (characterized by the condition that ri{X,j) = e) and the 
points for which ri{X,j) = e + j — A, and the second piece that consists of an infinite number 
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FIG. 7: A plot of the eigenvalues of the latticized Hamiltonian for 800 lattice sites with e = 0.001. 
Energies are plotted versus the variable kp = 27rp/800, for p = 0..799. 

of terms having essentially exactly the same form as the Hamiltonian for A = which has 
a spectrum which goes like k"^. Clearly the growth of a part of the spectrum which behaves 
like a free field with a very different speed of light from the rest of the theory represents 
a major change in the spectrum and thus constitutes a proof that this set of Hamiltonians 
really represent time dependent physics. 

In an attempt to achieve a better understanding of how these two different free field 
theories are linked together, we diagonalized the Hamiltonian in Eq. ( |104| ) for lattices which 
initially, in Schwarzschild coordinates, cover the region < r < 20Rs, where Rs stands for 
the Schwarzschild radius. Fig.|^ shows a plot of the spectrum for a lattice with 800 lattice 
sites, 5 = 0.025 and e = 0.001. The energy is plotted against the variable kp = |^ where 
p is an integer running from to 709 and represents what the lattice momentum would be 
for a very large number of lattice points. In this case we expect small corrections due to 
the finite size of the matrix corresponding to the kinetic term. We already noted that for a 
very large lattice we expect this plot to be proportional to k'^. In Fig.|| we plot this result 
scaled to have a maximum value of (27r)^ against k"^ over almost half the range of k in order 
to show that even for a finite lattice the agreement with our expectations is quite good. 

Next, in order to exhibit the difference between this situation and that for finite A, we plot 
the eigenvalue spectrum of the Hamiltonian for a lattice of 1200 sites, e = 0.001, 6 = 0.016 
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FIG. 8: We have scaled the energy to run from to 27r^/800 vs square of momentum for kp = 
27rp/800. The plot is limited to a bit less than half of the total range to show how small the 
deviations are from what we would expect in the infinite volume limit. 

and A = 300. This is a situation which corresponds to having about 300 lattice sites on 
the curve r = e and 900 points extending from r = e + 1 out to a distance of 20 times 
the Schwarzschild radius. Fig.^ shows two curves. The first is a plot of the lowest 300 
eigenvalues of the kinetic term divided by e^/'^ and plotted against a momentum variable 
kp = 27rp/300, where we have scaled the eigenvalues to make them show up on the plot. 
The choice of momentum variable is motivated by the fact that, as we have noted, the first 
300 terms in the Hamiltonian would have a spectrum almost proportional to (1 — cos(A;p)) 
if we dropped the term linking them to the next 900 terms. The second curve is a plot of 
the next 900 eigenvalues of the Hamiltonian versus a momentum kp = 27r/900. This should 
be quite similar to the curve shown in Fig.|^ and it is. We should emphasize that while for 
the case A = this spectrum starts at zero, in the case A = 300 it doesn't; admittedly that 
fact is not obvious in this plot. 



A cleaner picture of what is happening in the low energy region is shown in Fig.|TO|, where 
we have chosen to replot both curves. This time we haven't rescaled the first 300 eigenvalues 
and we have greatly expanded the vertical scale so as to see what is happening near zero 
energy. What one immediately sees is that the branch of the plot which we identified as 
belonging to the states localized on r = e is in fact massless, but the second branch, which 
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FIG. 9: The lower curve represents a plot of the first 300 eigenvalues of the Hamiltonian for 
A = 300. These correspond to the number of points which lie on r = e at this time. This curve is 
scaled to make it visible and it is plotted vs momenta kp = 27rp/300. The second curve represents 
the remaining eigenvalues, which naively correspond to the fields lying between r = e and r = 20R. 
This is plotted vs the variable kp = 27rp/900. 

we thought should look like the spectrum of the A = theory actually starts from a small 
gap which lies just above the final value of the first branch. How can this happen? The 
answer seems to be that the term which links the two distinct pieces of the Hamiltonian 
mixes what would have been degenerate levels and splits them to produce a continuously 
rising spectrum. In particular this would suggest that the lowest energy states we thought 
of as completely localized at r = e are in fact linear superpositions of such states and low 
energy states of the A = problem which are not localized at all. Clearly, it will require more 
work to convert these results into a better understanding of how and where the information 
which enters the horizon is really stored. 

Having shown that analysis of the spectrum of the time dependent Hamiltonian proves 
that the dynamics is time dependent, we wish to conclude with a few remarks about the 
scenario this picture suggests for the last moments of black hole evaporation. Clearly, the 
picture strongly implied by our discussion is that during the evaporation process a remnant 
is formed. This remnant represents almost zero-energy degrees of freedom which can be 
associated with the region r = e, since it is only by including these degrees of freedom that 
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FIG. 10: A plot of the two pieces of the eigenvalue spectrum of the A = 300 Hamiltonian on an 
expanded scale to show that the r = e and e < r < 20R mix and produce a gap. 



we preserve the unitary evolution of the system. The question therefore is what happens to 
these degrees of freedom when the black hole evaporates. 

One argument against a remnant would be that as holes evaporate these remnants should 
announce their presence in some dramatic fashion. Since our treatment of the problem 
encodes the changes in the metric into the changes in the coefficients appearing the field- 
theory Hamiltonian we can approach this question by asking what really happens to these 
coefficients during the final, rapid evaporation phase. If the coefficients of the Hamiltonian 
freeze into the form they had shortly before the point at which the methods used in this paper 
break down, then we see that finally the system will be composed of two very weakly coupled 
subsystems and the information stored in the remnant will not burst out, but at best dribble 
out over very long times. Obviously these remarks are not a proof of anything. At best 
they describe an alternative scenario for what might happen. What we wish to emphasize 
is that the Hamiltonian formulation of the problem coupled with lattice techniques gives 
one a new way of probing these issues. Hopefully a better understanding of whether all the 
information comes out will emerge as one studies this problem in greater detail. 
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XVII. SUMMARY 



In the preceding sections of this paper we argued that a consistent Hamiltonian formula- 
tion of the theory of a massless scalar field in the background of a Schwarzschild black hole 
exists, but at the expense of having an explicitly time-dependent Hamiltonian. We then 
reviewed the familiar discussion of the moving mirror problem as a way of emphasizing that 
a perfectly understood system with a time dependent Hamiltonian can evolve unitarily and 
still some observers will think that they are in a thermal bath of quanta. The purpose of this 
review was to emphasize that when one is dealing with a non-equilibrium system and one 
imposes the ideas of equilibrium thermodynamics one can come to surprising conclusions. 

Following the discussion of the moving mirror problem, we returned to the problem of 
the Schwarzschild black hole and used essentially the same techniques to show that, with 
reasonable assumptions about the initial state, our method leads to the usual result of 
Hawking radiation. 

There were two reasons for this discussion. First, we wanted to show that our treatment 
of the problem reproduces familiar results. Second, we used it to argue that formulating 
the theory on a spacehke shce at finite times clearly exhibits the fact that the Hawking 
radiation phenomenon emerges under rather general assumptions about the initial state on 
the quantization surface. We included a discussion of the case of a massless scalar field and 
an infalling refiecting mirror to sharpen this point. Assuming that the mirror is static for 
large times in the past, we can reasonably argue that it is proper to consider the system 
to be in the ground state of the theory, or in a state which differs from it by a finite 
energy. This naturally leads us to the usual late time Hawking radiation. Moreover, we 
saw that this Hawking radiation came from a place on the original surface of quantization 
which is far from the mirror and corresponds to the place from which null geodesies must 
depart so as to be scattered from the mirror within an exponentially small distance of the 
horizon. Clearly this is not a special point and its location depends upon when the mirror is 
allowed to move. Thus, we see that in order to get rid of the Hawking radiation by suitably 
modifying the initial state would amount to making a very strange ad- hoc assumption about 
the expectation values of field operators on the initial quantization surface. 

Having shown that our discussion leads to the usual picture of what external observers 
would measure we turned to the question of black hole entropy. Once again we discussed 
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the case of a black hole surrounded by a static mirror, but this time with the field theory 
restricted to the inside of the sphere. We then argued that from the point of view of an 
outside observer this system would look like a classical body with a temperature equal to 
the Hawking temperature and an equation of state which would imply, for an equilibrium 
system, an entropy equal to the Hawking-Bekenstein entropy Next, we looked inside the 
system and saw that it was never in equilibrium and that it was in fact always in a pure 
state. 

Finally, we identified the "information paradox" with the question of what is really hap- 
pening at the spacelike singularity r = 0. By discretizing our time-dependent Hamiltonian 
and studying the behavior of H{\) as a function of A we arrived at a very interesting picture 
which implied both that the modes at r = and the modes going out to large Schwarzschild 
r mix. This left us with the possibility that the endpoint of the evolution of the hole when 
it evaporates is a field theoretic system in which some of the information is stored in a 
very weakly coupled (but not decoupled) remnant and some has been squeezed out as it ap- 
proached r — 0. Clearly the questions raised in this section of the paper by far outnumber 
the results and these issues merit further study. 

In conclusion, we would like to reiterate the point made in the introduction to the paper. 
Our studies show that the theory of a massless scalar field in a black hole background is 
perfectly consistent with unitary time evolution and that somehow the theory resolves any 
supposed paradoxes in its own way. This led us to the conclusion that this semi-classical 
problem alone does not provide any smoking gun telling us what a correct theory of quantum 
gravity must look like. 
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